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Interferometers permitting the measurement of shape and altered distribution of fringes 
are usually used to compare either an unknown wave front with a known one or an unknown 
surface with a standard reference surface. Any error in the reference surface introduces a 
corresponding error in the results sought. This paper describes an absolute measuring 
interferometer and the associated mathematical operations necessary for analysis of the data. 
The reference surface is purely mathematical in nature and is therefore free from error. 
Deviations of converging wave fronts from perfect spheres, paraboloids or ellipsoids are 
readily measured without the use of tangible reference surfaces. 

The equation of the reference surface may contain one or more parameters whose value 
is sought. Thus, the eccentricity of the conicoid that best fits a mirror and deviations of the 
surface from the true mathematical curve are obtained simultaneously. The sensitivity of 
the test can be varied so that when testing large aspherical elements which depart very far 
from spheres the number and width of the fringes can be adjusted to any desired values. 



1. Introduction 

Many tests made by interferoiuetry consist in 
measuring the variation in separation of two wave 
fronts. Usually, the shape of one of the wave fronts 
is know^n so that the variation in the observed 
separation is a measure of tlie absolute shape of the 
unknown. The shape of the wave front is a function 
of the optical elements (lenses, mirrors, etc.) that 
produce it. Consequently, the aberrations of lenses, 
the shapes of mirror surfaces, the variations in densi- 
ties of fluids, etc., are obtainable from the measure- 
ment of wave-front shapes. 

Usually the known wave front is prochiced by a 
mirror of know sliape caUed the reference surface 
and tlie unknown wave front whicli is produced 
by the specimen is compared witli it. If tiie reference 
surface is imperfect, any error in it is reflected into, 
or transferred to, the measured sliape of the un- 
known. 

This paper describes a method of comparing wave 
fronts against a reference surface which is purely 
mathematical in nature. Since the Koesters double- 
image prism, which is the basic element of this 
interferometer, can be made with a high degree of 
perfection [1] ^ the results are practically free from 
instrumental error and the necessity i'or making 
reference mirrors and lenses of various foci is elimin- 
nated. The interferometer is quite stable, compact, 
and is easily operated by nontechnical personnel. 

This method of analyzing interferometer data, with 
the necessary modifications required to adapt it to 
the wave-front -shearing mterferometer ( = WSI), 
has been described in part 1 of this series of papers [2]. 



1 Figures in brackets indicate the literature references at the end of this paper. 



Without this development, neither the WST nor the 
wave front reversing interferometer ( = WKI), will 
yield unique results except when applied to wave 
fronts whicli are known (from some other type of 
test or source of information) to have revolution 
symmetry. 

2. Description of the Method 

To describe the principles used in this method of 
analysis its application to the testing of a simple lens 
will be used for illustration. The operations de- 
scribed here apply, without change, to the testing of 
lenses by autocollijuation (one conjugate at infinity), 
the testing of mirrors, and to all otlier tests where 
the results sought can be obtained froju a knowledge 
of wave-front shapes. This treatment is applicable 
to three dhnensional space but for simplification the 
mathematics will be limited to a plane. 

2.1. Symbols and Abbreviations 

The discussions and treatment of this problem 
require a considerable number of symbols. To 
provide a ready reference, all quantities used are 
defined below. 

a = diameter of the circular aperture of 
the lens or mirror under test. 

B = dividing plane of the double-image 
prism. 
C and C" = images of conic section. 
Ci, C2, and C2= arcs of reference curves representing 
known wave fronts; (i.e., arcs of 
circles if a corrected lens or spherical 
mirror is considered and arcs of 
a conic section if a comcoid is being 
considered) . 

D = position of small, monochromatic light 
source. 

£^= center of curvature of prism base. 

e = eccentricity of conic section. 
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T, r, Ti, T[, T2, T2 = 

2A and ?;= 

Fand V - 

W, TFi, 1^2, and W^- 



F= ideal image point (conjugate of D in 

fig. 1). 
/^j^ 7^2, F' = images of F {Fi and F2 are imaged by 
the outer plane surfaces of the prism 
and F' is imaged by the dividing 
plane B.) 
Fc and Fc = centers of curvature of conic section 
corresponding to points of extreme 
curvature. 
h and A; = coordinates of F^. 

/i and A; = components of the line from F to E, 
perpendicular and parallel to the 
dividing plane, respectively. 
iV,.ff= observed order of interference at P^. 
PvandP<T=any chosen pair of superimposed re- 
ference points. (These points are 
equally distant [optically] from the 
light source and, therefore, fall on 
the same wave front, Their posi- 
tions on the wave front are indicated 
to a sufficient approximation by 
spots on the lens surface, by beads 
on a thin wire or string, or by lines 
on a scale.) 
P= radius of curvature of a conic section; 
minimum for e"^ positive and maxi- 
mum for e^ negative. 
>S = separation of reference curves Ci and 
C2 (or C and C) along the line ot 
sight from E through Py. 
reference points on a conic section, 
coordinates of E. 
perihelion of conic section, 
respectively, the wave front before it 
is incident on the lens or mirror, 
the part of the transmitted wave 
front that lies above the dividing 
plane, the transmitted part that 
lies below the dividing plane and 
the image of W2 with respect to the 
dividing plane, 
wave-front-shearing interferometer, 
wave-front-reversing intertero meter, 
coordinates of T. 
angle between the dividing plane and 

the optical axis, 
angle between the dividing plane and 

the radius vector to py. 
separations between Ci and Wi along 

the radius vector to p^. 
separation between C2 and W2 along 
the radius vector to P^, (or separa- 
tion between C2 and W2 along the 
radius vector to Py). 
€ = angle between two images of conic 

section. 
€^ = angle between Ci and C2 at their 
point of intersection with the divid- 
ing plane. 
X = wavelength of light. 
ju = distance from the extension 01 the 
dividing plane to the axial center 
of the lens (fig. 1). It is positive 
when the dividing plane is above the 
axis of the lens and negative when 
below. 
V and (7 = distances from Py and P„, respectively, 
to the optical axis. These distances 
are measured along the arc of the 
reference curve, C1C2. They are 
positive when the corresponding 
points are above the axis of the lens 
and negative when below. 
P, p\ Pcu and pc2=radius vectors from E to T, T", Ti, 
and 7^2, respectively. 
py and Pa = radius vectors from E to P^ and P^. 
respectively. 
«^ = the sum of the squares of the devia- 
tions (defined by eq (13)). 
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2.2. Apparatus 

A simple interferometer for testing lenses, without 
a reference surface or element for comparison, is 
shown in figure 1. It consists of a small monochro- 
matic light source at D^ the lens to be tested and a 
double-image prism [1]. The lens received the 
diverging wave front, W (spherical when three dimen- 
sional space is considered and circular for two dimen- 
sional space) and changes it into a converging wave 
front. If the lens were perfectly corrected for the 
indicated position of the source, the converging wave 
front would be an arc of a circle (curve C1C2) and the 
light would form a perfect image (airy disk) of the 
source at the ideal image point F. The wave front 
produced by an imperfect lens (greatly exaggerated) 
is indicated by curve W1W2 and its deviation from 
the ideal wave front, C1C2, is a measure of the abbera- 
tions of the lens. The circle that represents the ideal 
wave front exists in mathematical form only and its 
parameters are adjusted statistically or otherwise 
chosen so that it fits the real wave front in accordance 
with some arbitrarily selected specification. Methods 
of choosing the reference circle will be described later. 

In practice the prism is mounted so as to pivot 
about E, and Ek placed near F. The semi-reflecting 
dividing plane, B, is adjusted so that its extension 
cuts the circular lens along a chord at a chosen dis- 
tance, iJL, from the axis of the lens. Each new adjust- 
ment of the prism (change in value of fx) introduces 
and additional arbitrary value for e^,. The quantity 
jLi is positive when the dividing plane is above the 
axis of the lens and negative when below. 

The effect produced by this prism on the wave- 
front is as follows: the wave front W, after division 
by wave-front division into two parts, W] above the 
dividing plane and W2 below, is further subdivided b3^ 
amplitude division at B into two beams each, making- 
four beams altogether, two arriving at each of the 
points Fi and F2. One of the beams that arrives at 
each of these points suffers two reflections and the 
other a single reflection. Consequently, a mirror image 
of everything seen below the dividing plane appears 
above it and vice versa. The arc C2 (fig. 2A) appears 
as C2' . The center of curvature of C2, [3] located at F, is 
imaged at F^ which is the center of curvature of 62^ 
The part W2 of the wave front below the dividing 
plane appears as W2', so that Wi and W2 appear 
inverted with respect to each other and superimposed 
in the overlapping region TI^2', where fringes of inter- 
ference (fig. 2B) are observed. 

The two halves of the prism are slightly rotated 
about an axis normal to the dividing plane (as de- 
scribed in references 1 and 6). A wave front having 
symmetry with respect to the dividing plane would 
therefore give rise to straight fringes normal to the 
line representing the dividing plane in the field of 
view. If the two halves of the prism are not ro- 
tated relative to each other the fringes are straight 
and parallel to the dividing plane. 

2.3. Equations of Observation 

The separation of wave fronts, Wi and TI'2' (fig. 
2 A), is the optical path difference pv—pn which is 
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Figure 1. Optical arrangement for testing a lens for a finite object distance. 
This arrangement requires the use of polarized light because of polarized effects in the prism. 



Figure 2. A. Illustration of image formation 
by the double-image prism. 
B. Interference fringes obtained 
with the assembly of figure 1. 




equal to the product XA^^^. The corres ponding sepa- 
ration of the reference circle C\ and G' at Py is, 
approximately [2]. 

S={u—ix)€^=€^{u—a). (1) 

The two quantities in parentheses are the distances, 
along the circle, from I\ and Pa, respectively, to the 
divi(Ung plane. Since P^, in figure 2, is below the 
axis of the lens, a is a negative quantity. 

It is apparent from figure 2A that the relation 
between the four separations, 8^, d^ S, and XA^^^ is 
given by the equation 

8.+S=8a+\N,.. (2) 

See section 2.1 for definitions of these quantities. 
On substituting for Sits value from eq (1) and trans- 
ferring terms we have 

8y—8a=\Npa—efji(v—fx) (3) 

whicli are the basic equations of observation, provid- 
ing relationships between the deviations of the wave 
front from the abstract curve C1C2. 

2.4. Distribution and Number of Reference Points 

If the surfaces of the lens are figures of revolution, 
the optic axis is well defined, and can be accurately 
located by observing the fringe pattern while chang- 
ing the value of jli. When the dividing plane coin- 
cides with the optic axis the fringes will be straight 
regardless of aberrations in the lens. If the lens 
surfaces are not figures of revolution the changes in 
curvature and spacing of the fringes will indicate tlie 
most probable position of the axis. 

The reference points are chosen in the aperture of 
the lens, along a straight line normal to the dividing 
plane of the prism and through the optical axis, 
regardless of whether the axis is or is not centered in 



the aperture of the lens. The positions of tliese 
points are measured from the optical axis. They 
are positive when above and negative when below 
the axis The points are equally spaced and their 
separation is defined as unity. An ochl or an even 
number of reference points, symmetrically placed 
above and below the axis, may be used An o(hl 
set, therefore^ would be located at integral units of 
length from the axis, as in figure 3-A, while an even 
set would be placed at odd integral half units of 
length from the axis, as shown in figure 3-B. If a 
point at the center is desired one must choose an odd 
number of points. Figure 3-C illustrates the ar- 
rangement of an odd set when the optic axis is not 
centered in the aperture. The interferometer is 
sufficiently sensitive to detect deviations of the po- 
sition of the axis from the center of the aperture in 
most lenses. 

Any number of points may be selected along the 
principal diameter (the diameter that is normal to 
the dividing plane) by adjusting the separation of 
adjacent points. If the mirror or lens to be tested 
is unsymmetrical and/or irregular the test should be 
made with more points than when a smooth wave 
front is formed. 

Even when the wave front is smooth and symmet- 
rical about the axis, many reference points will 
sometimes be required. The number of fringes in 
the field of view depends upon three quantities. 
These are: (a) The magnitude of the wedge that is 
built into the prism, (b) the angle betw^een the wave 
fronts at the dividing plane, and (c) the relative 
shapes of the two images, Wi and W2'? of the wave 
front. The effects of these three quantities are more 
fully described as follows: 

(a) The built-in wedge of the prism produces an 
invariant wedge between the wave-front images. Its 
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Figure 3. A, B, and C are, respectively, the distribution of the reference points when the num- 
ber is odd, when the number is even, and when the optic axis is not centered in the lens. 



vertex is perpendicular to the dividing plane. Since 
the component of this wedge, along the line of refer- 
ence points, is zero, its effect on the interference at 
the reference points is to add a constant to the 
observed orders of interference at all reference points. 
If the vertex of the built-in wedge is adjusted to 
coincide with the line along which the reference 
points are to be selected, this wedge does not affect 
the order of interference at the reference points. It 
does, however, ensure that the fringes never get too 
broad for precision readings. 

(b) The angle, e^, between the wave fronts at the 
dividing plane is equal to the angle between lines 
drawn from points F and F^ to the vertex of this 
angle (see fig. 13, A). This angle is adjustable by 
varying the position of irrelative to E, by displacing 
the prism laterally. Thus the component of fringe 
widths perpendicular to the dividing plane and the 
order of interference at an}^ chosen reference point 
can be adjusted to any desired value. If the two 
wave fronts, Wi and TF?', are identical in shape 
(i.e., if wave front W1W2 is spherical) the difference 
in order of interference between all pairs of adjacent 
points is approximately constant and adjustable at 
will. 

(c) If the wave front W1W2 is not spherical the 
two images of it, Wi and ^V, will not be identical 
unless 11/^1^2 has revolution S3^mmetry and the value 
of fjL is zero. When testing a wave front that de- 
parts very far from a sphere, as, for instance, when 
testing an //I parabola of large aperture, the mini- 
mum number of fringes that can be obtained by 
adjusting €^ depends upon fx. Since, in operation 
of the instrument, pairs of points are brought into 
coincidence for conditions other than fji=0, and since 
the absolute magnitude of fi depends upon the unit 
separation of points, it follows that the choice of the 
number of reference points determines the minimum 
attainable number of fringes. If the number of 
reference points is too small the minimum attainable 
number of fringes may be too large for accurate 
readings. Consequently, a large number of points 
may be required to obtain the optimum reading 
condition. Since the sensivity of the interferometer 
depends upon the ratio /d/a, the sensivity also depends 
upon the number of reference points. The variation 
of sensivity with the ratio /jLja and its effects on the 



required number of reference points are fully ex- 
plained under ^'Sensitivity'^ in part 1 of this series [2]. 

2.5. Solution of Wave-Front Deviations 

We shall select (fig. 4) quite arbitrarily for this 
discussion, seven reference points along the principal 
diameter of the lens and with their spacings sub- 
tending equal angles at E. These points will be 
considered as lying on a common wave front at the 
time it emerges from the lens, and their separations 
will be defined as one unit of length. Since the 
points are observed from position E they may be 
indicated with small beads on a fine wire that is 
stretched across the face of the lens (fig. 4A). 

If the dividing plane is adjusted to include one of 
the inner points (any point except an end one) as in 
figure 4C, some of the points will appear to coincide 
in pairs. If it is adjusted to bisect the angle be- 
tween any two of the inner points, as in figure 4B, 
some of the points again appear to coincide in pairs. 
The points will coincide in pairs when fx is an integral 
number of lialf units of length but will not coincide 
for other values of n. 

We may substitute, in eq (3), the observed values 
for Nj,a and (^— m) at each pair of points, successively. 
This yields as many equations of observation as 
there are pairs of coinciding points. We may then 
change the value of ju, thereby changing the com- 
bination of paired points, and obtain additional 
equations. Each new value for ijl (i.e., At=0,±0.5 
dzl.O, . . .) introduces several additional equations 
but only one additional unknown — i.e., a new value 
for €^. In this manner, we can obtain more equa- 
tions than we have unknowns; and it would seem 
that a statistical analysis of these equations would 
permit an evaluation of the several deviations (values 
for the 5^s). However, such is not the case. The 
relation between these equations are such that the 
resultant number of independent equations is al- 
ways less by three than the number of unknowns. 

When the wave front is completely unknown, the 
minimum number of observation equations that 
permit a unique solution of the 5's require two sets 
of fringes (two values for fi). Each adjustment of 
the interferometer introduces a new value for €„ 
because the value of e» can not be retained from 
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Figure 4. A. Illustration to show the manner of locating the reference points, Pj/, by means of 
heads on a fine wire. 
B. and C. Distribution of points and images of points for ix=—0.5 and fx=—l, 
respectively. 



one adjustment to another. For simplicity we will 
use the two adjustments (2 sets of fringes) that are 
obtained witli ju=0 and /x=0.5. These are (for the 
7 chosen reference points) for /x^O 

5i=5_i + XiVi._i-€o, (4) 

62=5_2+XiV2.-2-2eo, (5) 

53=5_3 + XiV3,-3-36o, (6) 

and for ^=0.5 

5i = 5o+XA^i.o-0.5€o.5, (7) 

52=5_i + XA^2.-i-1.5eo.5, (8) 

53=5_2 + XAr3._2-2.5eo.5. (9) 

Equations 4 through 9, form a set of 6 simultaneous 

equations having 9 unloiowns (7 5's and 2 e^t's); 

consequently, we either have an excess of 3 ujiknowns 

or a shortage of 3 eqiuitions for obtaining a unique 

solution. It is quite obvious that we cannot evaluate 

deviations from an unick^ntified curve. The persistent 

shortage of three equations is associated with the 

three conditions necessary to define an arbitrarily 

chosen 3 -parameter reference circle. 

Tlie method used in combining the three condition 
equations with the set of observation equations is 
fully described under the heading '^Specifications for 
the Reference CircJe" that appears in part 1[2]. 
Since the details of the application to the WRI are 
so similar to those given for the WSI, this treatment 
will be limited to details that are necessary for a 
clear understanding of the WRI. Since three points 
fully define a circle we may define the reference 
circle, C1C2, as that particular circle which passes 
through an}" three of the chosen points, by equating 
the corresponding 5^s to zero. Thus, if we impose 
the three conditions, 



^3 — ^0=5- 



= 



(10) 



we are requiring that the reference circle shall pass 
through points I\ Pq, and P_3, as indicated in figure 
5B. This may be considered as either a reduction in 
the number of unknowns by three or an increase of 
three in the number of independent simultaneous 
equations. In either case we get as many equations 
as unknowns so that a unique solution of all un- 
knowns is possible. Any three of the chosen reference 
points may have been selected to fall on the circle. 



Tlie simplest manner of specif3dng and computing 
the reference circle (called the '^Method of Coinci- 
dence^O may not produce a good fitting circle; but 
this does not render any less obvious the location of 
the high and low places on the wave front. Con- 
sequently, an optical worker would have no trouble 
locating tlie high and low places for fuilher polishing 
to reduce the deviations. 

A better fitting reference circle is obtained by ap- 
plying the method of averages [4]. Figure 50 shows 
the same wave front with a circle of reference that is 
determined by this method. The method of averages 
allows a variety of groupings of points and assign- 
ments of weights. For instance, since the areas of 
the wave fronts that are associated with the several 
points increases with v we might multiply aU 5's in 
eq (10) by the quantity (^+1) so as to assign unit 
weight to the axial point, 1\, and corresponding 
higher weights to other points, thus distributing the 
weight more nearly proportional to the areas of the 
several zones. 

A good fit, based on the method of averages, is 
obtained by imposing the following simple conditions : 



63 +^2 =0 

d, +do +5-1 = 
6_2+5-3=0. 



(11) 



The circle that fits the wave front shown in figure 5, in 
accordance with the conditions of eq (11) is shown in 
figure 5C. Equations (11), when added to the six 
equations of observation (eq (4) through (9)), again 
supply the necessary additional independent equa- 
tions for a complete solution of all 5's. 

The best fitting circle, shown in figure 5D, is de- 
fined bv a third set of three conditions, 



50 d0 d0 



where 



bdo deo deo.5 



= 



(12) 



0^ J2 di=8U + 8UW-r + dl+dl+5l+dl (13) 
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Figure 5. A. Illustration to show the distribution of the reference points relative to the prism. 
B, C, and D. The positions of the reference circles relative to the wave front when 
the reference circle is chosen in accordance vjith the method of coincidence^ 
the method of averages, and the method of least squares, respectively . 



which requires that the sum of the squares of the 
deviations be a minimum. Differentiation with 
respect to any other three of the nine unknowns 
(i.e., 7 5's and 2 e^^^s) will yield the same set of condi- 
tion equations. 

In order to determine the equations of condition 
corresponding to eq (12), it will be convenient 
to solve for the 5/s. These are : 

5o=5o (obvious) (14) 

5i = 5o+XA\,o-0.5eo.5 (from 7) (15) 

5_i=5o+X(A/'i.o-iVi.-i)-0.5eo.5+€o (from 4 and 15) 

(16) 

52=5o+X(M.o-M,-i + A^2.-i) 

-2.0eo.5+eo (from 8 and 16) (17) 

5_2 = 5o + X(A'l.O-M,-l+A^2.-l-A^2,-2) 

-2.0eo.5 + 3eo (from 5 and 17) (18) 

53 = 5o + X(iVi.O-M.-l+A\-1-iV2,-2+iV3.-2) 

-4.5€o.5+3€o (from 9 and 18) (19) 

6_3 = 5o + X(A^l,0-M.-l+iV2,-l-iV2,-2 + iV3,-2-iV3,-3) 

-4.5eo.5 + 6€o (from 6 and 19) (20) 

On performing the differentiations indicated in eq 
(12), and substituting the corresponding values for 
these differentials from eqs (14) through (20), we 



obtain the set of condition equations for the method 
of least squares. These are: 



53+52+5i+5o+5_i+5_2+5_3-0 
6_i+52 + 3(5_2+53)+65_3-0. 

(5i + 6_:)+4(52+5_2) +9(63+5,3) = J 



(21) 



A.ny one of the sets of conditional equations 
(given by eqs (10, (11), or (21) combined with the 
observation equations (4) through (9)) permits a 
complete solution of all unknowns. The operations 
that lead to a solution of these equations are per- 
formed as follows: (1) B3' adjusting the interferom- 
eter so that /z=0 (i.e., with the dividing plane 
cutting the wave front along a line through the 
axis of the lens; (2) substituting the observed 
values for (^— m) and the observed orders of inter- 
ference at each of the three pairs of coinciding 
points into eq (3) thus obtaining eqs (4), (5), and 
(6) : (3) readjusting ^x to equal 0.5 (fig. 4B) : (4) 
again substituting the observed quantities Nj,a, v 
and \x into eq (3), thus obtaining eqs (7), (8), and 
(9) : (5) solving for the 9 unknowns in the 9 linear, 
simultaneous equations (eqs (4) through (9) and any 
one of (10), (11), or (21)). 

2.6. Symmetrical Systems 

Lenses and mirrors that have been polished 
mechanically produce symmetrical wave fronts 
when tested on axis. For such symmetrical elements 
all 5's on one side of the axis of symmetry are equal 
to corresponding 5^s on the other side and for ^t=0 
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{he fViii^(*s of inl(M'IVroii('(> will be straight. All 6's 
with negative subscripts can, therefore, be replaced 
with the corresponding positive subscripts. In the 
case of seven reference points, chosen here for illus- 
t ration, the left hand side of eqs (4), (5), and (6) 
becomes zero and the right liand side of each of 
these equations represents a solution of eo. Since 
€n if^ superfluous, the adjustment of /x to zero for a 
symmetrical wave front serves only as a check on 
symnu^tiy \ Consequently, if symmetry is observed, 
ecjs (4), (5), and (6) are exempted from the group 
of simultaneous equations. The elimination of 
negative subscripts reduces the number of unknowns 
so that a single adjustment (/x=0.5, eqs (7), (8) 
and (9) of the interferometer is sufficient for a 
complete solution. Equations (4), (5), and (6) 
are ehminated and eqs (10, (11), and (12) for the 
three methods, are reduced, respectively, to 8^ = 



^0=0, d-s + d2 = 28i+8,=0, and 



d(/) _5</> 

^60.5 d5o 



=0- 



The min)l)er of e({uatioiis in the complete set is 
reduced from 9 to 5 with a corresponding reduction 
in the mimhei' of unknowns. 
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2.7. Illustration 

The process described above is ilhistratcMl 
following application to the testing of a 12-in. aper- 
ture, 15-ft focal length asti'onomical refractor. 
Figure 6 shows the optical arrangement for testing 
with one conjugate at infinity. Figure 74 is a 
photograph of the fringes taken with tlu^ dividing 
plane coinciding with the optic axis of the lens 
(^1=0) and 7B shows the fringes when At=/2 unit 
( = 2 cm). Eight reference points were used for this 
test. These were located at distances of ±2, ±6, 
±10, and ±14 cm from the optical axis. The unit 
separation is 4 cm. Therefore the absohite vahies 
for p and a are 0.5, 1.5, 2.5, and 3.5 units. Table 1 
shows the observed orchu's of interference for the two 



sets of fringes at the pairs of rcd'erence points indi- 
cated by the v and a columns. The straight ness of the 
fringes in figure 7A shows the wave to be symmetrical 
about the dividing plane. When straight fi'inges 
are obtained for several azimuthal positions of the 
lens (by rotating it relative to the prism about its 
optical axis) the test indicates the lens to be a figure 
of revolution. This means that 8^=8^ and eqs (4), 
(5), and (6) are not required. The solution for the 
5's can, therefore, be obtained by using only three 
observation equations and two equations of condi- 
tion. The equations of condition may be based on 
either one of the three criteria: i.e., method of 
averages, method of coincidence, or method of 
least squares. 

We shall arbitrarily choose to use the method of 
averages and select the reference circle bv requiring 
that 



00.5 + ^1.5=0 
^2.5+53.5 = 0- 



(22) 



Tlie equations of observation are obtained by 
substituting corresponding values for v, or, /z, and 
A^^^, shown in table 1, in the columns under (ju=0.5), 
into (Hi (3). These are: 



5i.5=5o. 5+0.55 —1 eo.o 

52.5=5i. 5 + 0.20 -2 €0.5 
^3. 5=^2. 5— 0.35 —3 €0.5.^ 



(23) 



Equations (22) and (23) contain 5 unknowns, 4 5's 
and €0.5. The simultaneous solution of these equa- 
tions for the 5's yields the deviations of the wave 
front from a sphere. These values are —0.23, 
±0.23, +0.28, and —0.28, respectively, for 60.5, 
5i.5, 52.5, i^nd 53.5. Since the light passes through the 
lens twice in this test these deviations are twice 
the values that would have been obtained for one 



OPTICAL FLAT 




FicuRE 6. Optical arrangement for testing a lens with one 
conjugate at infinity. 
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Figure 7. Photographs of fringes with 12- 
inch aperture, f/15 astronomical objective 
when (A) ju=0, and (B) ix = 0.5. 



Table 1. Data from figure 7 for 
analysis of a 12-inch aperture 
astronomical objective 



AX=0 


M=0.5 


V 


(T 


N,a 


V 


cr 


iV.a 


0.5 
1.5 
2.5 
3.5 


-0.5 
-1.5 
-2.5 
-3.5 


-0.10 
-0. 45 
-0.70 
-1.00 


0.5 
1.5 
2.5 
3.5 


0.5 
-0.5 
-1.5 
-2.5 


0.00 

0.55 

0.20 

-0.35 



transmission through the lens — as for hght from 
a celestial star. The deviations for one transmis- 
sion are represented graphically in figure 8. The 
cnrve indicates that when this objective is used 
for astronomical work the wave front departs from 
a perfect sphere by approximately one-third of a 
wavelength. The light used for this test was the 
green line of the mercury spectrum (X= 0.546 
microns). 

If eq (3) is used in the testing of a parabolic mirror 
the computed deviations will be the deviations of the 
mirror from a perfect sphere. Usually, when testing 
a parabolic mirror one wishes to know how much the 
mirror differs from a perfect paraboloid. If the 
mirror is ellipsoid a 1 one wishes to know both the de- 
viations and the eccentricity. The eccentricit}^, 6, 
enters as an additional parameter (eq (30) of ap- 
pendix A). 

The simultaneous evaluation of the deviations and 
the eccentricity is obtained with the method of coin- 
cidence by requiring the curve to pass through 4 
points instead of 3 as was the case (eq (10)) for a 
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Figure 8. Graphical representation of wave-front-shape com- 
puted frorn the photographs of figure 7. 

3-parameter reference curve. If the method of aver- 
ages is to be used the reference points are grouped 
into 4 groups yielding 4 equations of condition in- 
stead of the 3 shown in eq (11). Similarly, with the 
method of least squares an additional equation 
(d0/d6=O) is to be added to the three equations under 
eq (12). If the surface is known, or can be assumed, 
to be a figure of revolution about its axis, the number 
of condition equations is reduced by 1 , as for the cir- 
cular reference curve. 

Figure 9 is a photograph of fringes taken with a 
12-inch aperture, 84-inch focal length, mirror that 
was claimed to be a nearly perfect paraboloid. Fig- 
ure 10 shows the optical arrangement used for making 
this photograph. The computed shape of this mir- 
ror indicated that it conformed more closely to an 
ellipsoid than to a paraboloid. The general equation 
for an ellipse (described in appendix A, eq (31)) was 
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Figure 9. Photograph of fringes taken with a 12-inch aperture 
fl7 parabolic mirror. 




-LIGHT SOURCE 



Figure 10. Optical arrangement for testing co7icave mirrors at 
their centers of curvature. 




Figure 11. Photograph of two Koesters double-image prisms^ 
one mounted in a cell and the other untrimmed. 



llioii tried urid the computed value for the eceeii- 
tricity corresponded to an elKpse whose foci are sepa- 
rated by 89 meters. It was suhseciuently learned 
that tlie final test on this ''Nearly Perfect Parabolic 
Mirror" was made with a Foucoult knife edge test, 
with the source located 90 meters (assumed to be 
practically at infinity) from the mirror and the knife 
edge at the point of convergence or image of the 
source. This test, therefore, confirms the description 
of the mirror as being a "nearly perfect parabola" 
but even more nearly a perfect ellipsoid. 

Figure 11 shows photograplis of the prism, before 
and after it is fitted into a cell. This prism has been 
made in sizes from approximately 6 inches down to a 
fraction of an inch. The one shown here has a one 
inch aperture and is considered a most favorable size 
for many applications. 

2.8. Summary 

The interferometer described here for testing wave 
fronts is believed to be easier to operate than any 
previous interferometer yet developed for this pur- 
pose. Its operation is simple. Unskilled personnel 
can be trained to operate it after a short period of 
training. . 

This interferometer yields results that are abso- 
lute. It does not require a standard of reference. 
The absolute evaluation of wave-front shapes, how- 
ever, depends upon a unique method of analyzing 
the data. The analysis involves the development of 
a set of linear simultaneous equations. This requires 
some technical background knowledge. After the 
equations for a given assembly are developed, un- 
skilled personnel can compute wave-front shapes by 
substituting observed orders of interference into the 
equations and solving for the deviations from the 
chosen reference curve. The chosen reference curve 
will usually be a circle but more complicated curves 
may be used. 
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FiGiire 12. A. The family of curves representing the conic section for various indicated values of e^ (e = eccentricity), 
B. Illustration for describing the transformation from rectangular to polar coordinates. 



3. Appendix A. Optical Path Difference 
Between Any Two Points on a Conic Sec- 
tion to Any Point in the Neighborhood of 
the Center of Curvature of the Section at 
Its Vertex 

The equation for the conic section in rectangular 
coordinates, with the origin at Fc, (a distance B to the 
left of the curve, see fig. 12A), and the x-axis coincid- 
ing with the axis of the conic, is given by the 2 param- 
eter equation 



y'+{l-e^){x-Ry+2E{x-B) = ^. 



(24) 



The family of curves, shown in figure 12 A repre- 
sents this equation for the various indicated values 
of e^. The parameter B is the radius of curvature 
of the curves at the vertex V (except for 6^=zb 0°), 
and e is the eccentricity. 

If we shift to polar coordinates (p, /5) with the polar 
axis making angle a with the axis of the conic (see 



fig. 12B) and the pole located at point E, which is 
displaced from Fc by a small amount u parallel to 
the polar axis and a small distance v perpendicular to 
it, we may obtain the polar equation for the conic 
by means of the transformation equations. 



y=u sin a-\-v cos a+p sin (a:+/3) 
X=^U(^OSa—V sin a + p cos (a:+/3). 



(25) 



The quantities u and v will never exceed the separa- 
tion of points E and Fc, which is always quite small 
relative to p, because these two points can always be 
superimposed to a high degree of accuracy. 

On eliminating x and y from eqs (24) and (25) : 
solving for p; applying binomial expansions to elimi- 
nate radicals and negative power terms; expanding 
sin {a-^fi) and cos {a-\-p) in powers of (a:+/5) [5]; 
and dropping all power terms in u and v higher than 
unity (because these terms are always insignificantly 
small), we obtain 



where 



J, (a,^) = E+l/S Re' {a+py-1/^SBe' (1+3^ (a+/3)«+. . 
^2 (a, ^)=C0Sp + l/2 e' ^ (a+^y + . . . 
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(26) 





Figure 13. A and B. Illustration of the relative positions of two images of a wave front with the WRI and WSI^ respectively. 



The term /i(a, (3) is the value of p when the pole is 
located at Fc. When testing optics of large aperture- 
to-focal-lengtli ratios the quantity j\{a, (3) must 
include higher order terms of the series than for 
smaller /-values. It will be found that, in some 
cases, the j\{a, /3) -series converges too slowly for 
practical application and that the rigorous, unex- 
panded value 



/i(«,^) 



_^[Vi- 



-e'^ s'm'^{a+ 13)- e^ COS (a+p)] 
[l-e'cos'(a+/3)] 



(27) 



is more practical. For this treatment we shall limit 
the angles to values comparable to those obtained 
with an//8 system. This limitation makes the last 
terms, indicated in each of the functions /i(a, (3) 
/2(a, 13) and/3(af, (3) of eq (26) negligible so that eq 
(26) becomes. 

p=R+l/SRe\a+p)^-u cos ^-v sin /?. (28) 
3.1. Equations for the WRI 

When the conic is used as the reference surface in 
a w-ave-front-reversing interferometer, the optical 
path difference, S (fig. 13A), is the difference between 



two radius vectors, pi^ and p2c, forming angles with 
the dividing plane, w^hich are ecj[ual in magnitude but 
opposite in sign. 

The two curves, C and C^ shown in figure 13A 
are images of the conic section. These two images 
are inverted by the interferometer, with respect to 
each other, about the dividing plane which is parallel 
to the polar axis. The two points 7\ and T/ are 
images of a common point on the original curve. 
Similarly T2 and T2' are images of a common point; 
but Ti and T? are images of different points that are 
at equal angular distances from the dividing plane. 
The radii vectors to Ti and 7\ are pic and p2c, respec- 
tively. 

If we define the perpendicular displacement of Fc 
from the dividing plane as —h, tlie corresponding 
displacement for Fc is +h. The component dis- 
placements of E from F^ as defined for figure 12, is 
u and V. The corresponding displacements of E 
from Fc is u and (v—2h). The value of pic is, there- 
fore, obtained by replacing p in eq (28) by pu. Simi- 
larly, since the sign a is reversed in the image, the 
value for p2c is obtained by replacing p, a, and v by 
P2c, —a, and {v—2K)y respectively. On making these 
substitutions, we obtain 

P2c-Pic^S=2h sin p-Re'a^{a'+^'). (29) 
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It is convenient to use the measured distances 
(p—fji), and Id instead of the angles ^ and a. Also, 
from figure 13A, it is obvious that 2h=eE, since h is 
quite small relative to R. On substituting for 2h, a, 
/3 (=sin /3, approximately), their respective equals, 
eR, iJiR~\ and (v—jjl)R~\ the value of S becomes 

S^6(.-M)-r7r-V(^-M)[M^+(^-M)-1 (30) 

which becomes for a paraboloid (e = l), 

>S=6(.-M)-/rV(^-M)[M^+(^-M)-1 (31) 

and for a sphere (^ = 0), 



S=e(v — fM). 



(32) 



Equation (32) is the same as eq (1). The use of 
a conic section instead of a circle for the reference 
curve, in the derivation of eq (3), leads to the intro- 
duction of the additional unknown, e, into eq (3) so 
that an additional condition equation becomes neces- 
sary for a unique solution, thus permitting the 
incidental evaluation of e. Since four conditions are 
required to define the chosen conicoid, four condition 
equations are needed. These may be: The assign- 
ment of arbitrary values (usually zeros) to any four 
of the 8^s (if the method of coincidence is used) ; the 
grouping of the reference points into four groups and 
requiring the algebraic sum of the 8^s in each group 
to be zero (if the method of averages is used) ; or, 

the addition of the condition equation, ^=0, to eq 

(12) (if the method of least squares is used). 

3.2. Application to the WSI 

When the conic is used as a reference surface in a 
wave-front-shearing interferometer, the optical path 
difference, S, is (p^—p), where p and p' are the radius 
vectors from the point E to points T and T^ (see fig. 
13B) on curves C and C\ respectively. Points T 
and T^ are images of a common point on the conic 
section. Curves C and C^ are sheared images of the 
conic section. The angle of shear is —a for C and 
a for C\ Therefore, the polar axis makes angles 
a and ~a with the axis of C and C\ respectively. 

The component displacements of E from Fc is the 
same in figure 135 as in figure 125. The value for 
p is, therefore, given by eq (28). The component 
displacements of £ from F^ is (u—2k) and (v—2h). 
Consequently, the value for p' is obtained from eq 
(28) by replacing p, a, u, and v by p, —a, (u—2k) 
and (v—2h), respectively. On making these substi- 
tutions and taking differences we obtain 

p^-p=S=2h sin p+2k cos I3-Re'a^(a'+I3'), (33) 



It will be noted that k does not appear in eq (29) 
(relating to WRI) but does appear in (33) (relating 
to WSI). The term 2kcos^ is approximately con- 
stant when 13 is limited to small values (as when 
testing an //8 system) but becomes significant for 
large values of 13 (as when testing an //I system), 
unless k is sufficiently small. If the WSI is adjusted 
to place the zero order of interference at the center, 
the magnitude of k is insignificant and the values of 
S for the two interferometers are identical. On sub- 
stituting for 2h, a, and /3 their respective equals, as 
given in section 3.1, above (cos 13=1, approximately), 
eq (33) becomes 

S=e(v-f.)+2k-e'R-'fji{p-lJi)[fd'+{v-}xy-] (34) 
which becomes for a paraboloid (^^ = 1), 

S=e(v-fji)+2k-R-'ui(v-fji)[lJL'+{v-uiy] (35) 

and for a sphere (e = 0), 

S=e(v-iJi)+2k, (36) 

It is inconvenient and unnecessary to reduce k to 
an insignificant magnitude. It can be combined 
with other constants, as described in reference [2]. 
In fact, the constant ^ is a part of the constant r of 
reference [2]. Unfortunately, the sign convention 
for measuring p and v is different in this paper (part 
2) from that of part 1 [2]. The constant, r, com- 
puted from eq (3) of reference [2] in combination 
with eq (36) above, becomes r=Qr + eyL—2k. Since r 
is a constant and serves as a parameter only, its 
form does not affect the problem. 
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Fractional factorial designs for experiments with factors at 
two and three levels, W. S. Connor and S. Young, NBS 
Applied Math. Series 58 (1961) 40 cents. 

This publication contains fractional factorial designs for use 
in experiments which investigate m factors at 2 levels and 
n factors at three levels. The grand mean, all main effects, 
and all two-factor interaction effects are estimated. All 
higher-order interactions are assumed negligible and their 
absence allows estimation of the error variance. A design 
has been constructed for each of the 39 pairs {m,n) included 
from m + n = 5 through m-h?i=10, (m,n?^0). 

An evaluation of Kacser's second order Born approximation 
to the bremsstrahlung differential cross section, i\. 8. Ofell, 
NBS Tech. Note SI {PB161582) {1961) 75 cents. 

The second order term, as derived by C. Kacser [Proc. Hoy. 
Soc. A253 (1959)], of the Born approximation series for the 
bremsstrahlung cross section differential with respect to 
photon energy, i)hoton direction, and final electron direction 
averaged over initial and summed over final polarization 
states has been put in a form that admits to numerical 
evaluation for the coplanar case. The results are valid for 
relativistic as well as non-relativistic incident electrons. 
Completely general (noncoplanar) expressions are included for 
the first order (Bethe-Heitler) and second order terms. How- 
ever, the latter has not been given in a form admissible to 
numerical evaluation. 

Tabulated values are given for the case of the incident 
electron kinetic energ3^ = 500 kev, photon energy = 450 kev, 
and photon direction 20° from the incident electron direction. 
For this case the second order term increases the differential 
cross section by about 50 percent for Al and about 200 percent 
for Au for the various angles of electron emission. A rough 
approximation to the integrated (over final electron direction) 
cross section is included and a comparison is made with 
experimental values. The inclusion of the second order term 
brings the theoretical cross section closer to the experimental 
values. 

An approximate full wave solution for low frequency electro- 
magnetic waves in an unbounded magneto-ionic medium, 

W. C. Hoffman, J. Research NBS 66D {Radio Prop.) No. 1, 
107 {J an. -Feb. 1962), 70 cents. 

Maxwell's equations in an anisotropic inhomogeneous medium 
are transformed by means of the Stratton-Chu formula into 
a vector integral equation which couples the various electric 
field components. In case the hypotheses of far-zone field 
and low-frequency electromagnetic waves apply, this vector 
integral equation can be approximated by a system of un- 
coupled scalar integral equations. This imphes an approxi- 
mate equivalence between the original vector integral equa- 
tion and a system of modified scalar inhomogeneous Helm- 
holtz equations. The conditions under which the system of 
uncoupled scalar integral equations can be solved by Neu- 
mann series are discussed, and the first three terms of the 
Neumann series are given explicitly. 

Current on and input impedence of a cylindrical antenna, 

Y. M. Chen and J. B. Keller, /. Research NBS 66D {Radio 
Prop.) No. 1, 15 (Jan.-Feb. 1962) 70 cents. 

The current on a finite antenna is expressed as the sum of a 
current emanating from the gap and two currents reflected 
from the ends. These currents are determined for a perfectly 
conducting hollow pipe of circular cross section. It is excited 
by an electric field parallel to the axis applied across a gap 



of finite width which encircles the antennas. They are also 
determined for a thin antenna of any cross section. From 
the results the current on and input admittance of the corre- 
sponding finite antennas are determined. It is shown that 
the thin antenna theory yields an incorrect result for the 
admittance because it ignores a boundar}^ layer effect near 
the gap. 

On the propagation of VLF and ELF radio waves when the 
ionosphere is not sharply bounded, J. 11. Wait, /. Research 
NBS^^H {Radio Prop.) No. 1, 53 {Jan.-Feb. 1962) 70 cents. 
Employing an idea of Brekhovskikh, an expression for the 
reflection coefficient of a continuously stratified ionized 
medium is derived. The result is in the form of a series 
whose first term is a Fresnel-type coefficient and succeeding 
terms account for the finite thickness of the transition layer. 
This result is then fitted into previously developed theory for 
propagation between a spherical earth and a concentric 
ionosphere. 

On the theory of wave propagation through a concentrically 
stratified troposphere with a smooth profile. Part II. Expan- 
sions of the rigorous solution, 11. Bremmer, J. Research NBS 
66D {Radio Prop.) No. 1, 31 {Jan.-Feb. 1962) 70 cents. 
This part concerns the height-gain differential equation in 
order to obtain a series for the complete solution which starts 
with the extended W. K. B. approximation discussed in part I. 
The coefficients of this equation depend for each mode 
amongst other things on the parameters Aj fixing the refrac- 
tive-index profile. However, the explicit dependence on 
these parameters can only be given in terms of expansions 
with respect to {Ka)--'^ (A;oa = circumference of the earth 
divided by the wavelength). In turn these expansions are 
derived with the aid of other ones for the complex turning 
point connected with the height-gain differential equation. 
The final expansion for the solution of the differential equa- 
tion is substituted in the boundary condition at the earth's 
surface. This leads to corresponding expansions, with respect 
to (/coa)-2/3 of the quantity u X (a) and next of the eigenvalues 
I themselves. 

Reflection and transmission of radio waves at a continuously 
stratified plasma with arbitrary magnetic induction, J. R. 

Johler and J. D. Harper, Jr., J. Research NBS 66D {Radio 
Prop.) No. 1, 81 {Jan.-Feb. 1962) 70 cents. 
A flexible theoretical model plasma which can be deformed to 
fit most measured electron-ion-altitude profiles is employed 
together with available geophysical data on the ionosphere to 
evaluate reflections and transmissions during quiescent and 
disturbed propagation conditions. The reflections and trans- 
missions in the ionosphere are determined rigorously with the 
aid of the classical magneto-ionic theory. The complex 
indexes of refraction of the medium are deduced, and a 
couphng in the plasma between ordinary and extraordinary, 
upgoing and downgoing modes of propagation is investigated. 
The corresponding reflection and transmission coefficients are 
then calculated, and certain phenomena which can be ex- 
pected as a result of the action of a solar disturbance on the 
reflection process are predicted. 

The disturbance of solar origin, investigated as an application 
of developed techniques, influences the reflected and trans- 
mitted LF waves in the lower ionosphere in a complicated 
manner. However, the high absorption phenomena exhibited 
by high frequencies does not seem to exist for the plasma 
profiles investigated with the classical magneto-ionic theory. 
The electron collision frequencies of the classical magnetic- 
ionic theory are modified to introduce a collision frequency 
proportional to the electron energy, and the changes necessary 
in the formulation of the classical theory as a result of such a 
consideration are presented. 
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On the diffraction of spherical radio waves by a finitely 
conducting spherical earth, L. C. Walters and J. R. Johler, 
J. Research NBS 66D (Radio Prop.) No. 1, 101 (Jan.-Feb. 
1962) 70 cents. 

The theory for the diffraction of spherical electromagnetic 
waves by a finitely conducting spherical earth was developed 
from Maxwell's equations by Watson [1918] and the intricate 
computation details were later worked out by van der Pol 
and Bremmer [1936] as the now classical series of residues. 
Two aspects of this computation present considerable 
difficulty, especially at low^ frequencies: (1) the calculation 
of the height-gain factor which takes account of an elevated 
transmitter and/or receiver, (2) the evaluation of the special 
roots T^=Ts of Riccati's differential equation, (18/ dr 
— 252r+l = 0, near the circle of convergence, | S^r | =%. 
These analytic difficulties are avoided with the aid of m(^dern 
analysis techniques applied to a large scale electronic com- 
puter. Hankel functions of the first and second kind of 
order one-third and two-thirds are calculated by numerical 
integral methods and then used with iteration to solve 
Riccati's differential equation. The amplitude and phase of 
the spherical radio wave diffracted in the vicinity of the earth 
with various altitudes above the surface of the earth of both 
the transmitter and the receiver, are then calculated by a 
summation of the series of residues. 

The electromagnetic fields of horizontal dipole in the 
presence of a conducting half space, J. R. Wait, Cari. J. Phys. 
39, Xo. 7, 1017-1028 (July 1961). 

The problem considered is a horizontal electric dipole which is 
located above or below the plane surface of a conducting 
half-space. Expressions for the fields are obtained using 
three different approaches. The formulas developed are 
((uite simple and, taken together, the whole range of distances 
from the far-zone to near-zone is adequately covered. 

Simple calculus for all-dielectric interference filters, K. D. 

Mielenz, /. Opt. Soc. Am. 50, No. 10, 1014 (Oct. 1960). 
Chebychev polynomials are introduced into the matrix theory 
of first order all-dielectric interference filters consisting of 2rn 
or (2m-\-l) alternating high and low index films of equal 
optical thickness /? on each side of the low index spacer of 
thickness 2(3. For the amplitude and energy transmitted by 
such filters, simple closed expressions are derived that will 
render numerical results with greater facility than any other 
calculus known. 

Bounds for determinants with positive diagonals, E. V, 

Haynsworth, Trans. Am. Math. Soc. 96, 395-399 (1960). 
A new bound is given for determinants with ''dominant" 

11 
diagonals, that is those for which aii>Xj ^iJ (^=^1^ • - • > ^0, 

3=1 

and another is given for determinants whose diagonals are 
positive and satisfy the hypotheses of a theorem by A. 
Ostroski. 

The classical field theories, C. Truesdell and R. Toupin, 

Handb. der Physik, Springer, Berlin, III, No. 1, 226-793 

(1960). 

This paper presents an organized account of everything known 

concerning general continuous media. 

Resonance characteristics of a corrugated cylinder excited by a 
magnetic dipole, J. R. Wait and A. M. Conda, IRE Trans. 
Ant. Prop. AP-9, 330-338 (July 1961). 

Radiation from an axial magnetic current element in the 
presence of a corrugated cylinder is considered. It is indicated 
that the power radiated in a given mode depends on the 
surface reactance, the circumference of the cylinder and 
elevation angle. For certain values of the parameters 
particular modes are strongly excited, corresponding to a 
resonance condition of the circumferential (or spiral) surface 
waves. 

The calculation of energy levels of nearly symmetric rotors, 

S. C. Wait, Jr., J.Mol. Spectroscopy 6, 276 (1960). 
A recursion relation is developed for expressing the inverse 
of a power series by another power series. The derived 
relation is exact for any desired degree of approximation. 



An example of this recusion formula is given for use with an 
extended perturbation treatment of the nearly symmetric 
rotor. 

Theory of vibrational relaxation in liquids, R. Zwanzig, 

/. Chem. Phys. 34, 1931-1936 (1961). 

A new formulation of the theory of vibrational relaxation, 
based on Zener's semiclassical approximation, is presented 
here. The relaxation rate is shown to be proportional to the 
spectral density of the force exerted on the oscillator by its 
environment. The isolated binary collision theory is derived, 
but only with the condition that the collision frequency is 
much smaller than the oscillator frequency. This requirement 
is not satisfied in a liquid; we conclude that Litovitz's applica- 
tion of the isolated binary collision theory to liquids is not 
justified. A possible relation between vibrational relaxation 
and the self-diffusion coefficient in a liquid is discussed. 

Density fluctuations and heat conduction in a pure liquid, 

R. E. Nettleton, Phys. Fluids 4, No. 1, 74 (Jan. 1961). 
Two components of the heat flux in a pure liquid are identi- 
fied — one carried by sound waves and the other by diffusing 
molecules. Coupled relaxation equations are obtained for 
these components, both depending on the density gradient 
Vp, and this dependence is interpreted thermodynamically 
as a further coupling of these equations with the equation 
determining (d/dO(Vp). The latter can be immediately 
written down with the aid of Onsager's theorem and an 
inertial principle developed in a previous paper. The 
thermodynamic interpretation of these rate equations also 
leads to an explicit expression for the coefficient of (Vp)^ in a 
Taylor expansion, about equilibrium, of the Helmholtz free 
energy. This result is compared with similar free energy 
terms assumed to exist by other authors. 

Non-reciprocity and time-reversal in microwave circuits, 

H. A. Fowler, Academic royale de Belgique (Classe Des Sciences) 
V-46, No. 11, 963-977 (Dec. 1960). 

The relation between reciprocity of the microwave scattering 
matrix and a ''local" time-reversal symmetry of the scattering 
process is discussed, under the unitary condition. The 
theory of reciprocal junctions is extended to include non- 
reciprocal ferrite devices: the 7-circulator of Chait and Curry 
is analyzed in accordance with the point-group method. 

Tchebycheff approximations by ab^-fc, J. R. Rice, J. Soc. 
Ind. Appl Math. 8, No. 4, 691-702, (1960). 
In this paper the theory of Tchebycheff approximations by 
ab^-\-c and a6^ cos (do^dx) +c is studied where a,b,c, do and 
are parameters to be studied. A complete development of the 
theory is made for ab^' + c. The topics considered are exist- 
ence theorem, characterization of best approximations, 
approximations on subsets, uniqueness theorem, approxima- 
tion on finite point sets, limiting relations between approxi- 
mations on [0,1] and finite point sets and computational 
procedures for finite point sets. A really surprising phenom- 
enon of this theory is that best approximations do not neces- 
sarily exist on finite subset of [0,1], although best approxima- 
tions on [0,1] exist for every continuous function. The 
results for ab^ cos (do-\-6x) +c are very incomplete and only 
the rudiments of a theory are established. This is due mainly 
to the intractibility of the transcendental equations arising 
in the analysis. On the other hand some useful formulas for 
apphcations are obtained. 

Repulsion of energy levels in complex atomic spectra, R. E. 

Trees, Phys. Rev. 123, No. 4, 1293-1300 (Aug. 15, 1961). 
Rosenzweig and Porter have shown a "repulsion of energy 
levels" in spacing distributions determined from energy 
levels in complex atomic spectra. The present paper extends 
their work by showing that these spacing distributions can be 
determined from calculated positions of the levels in these 
spectra. Since calculated data are available for spectra 
where the observed data are scarce or incomplete, this par- 
tially overcomes limitations imposed by statistical inaccuracy 
when direct use is made of the observed data. The equiva- 
lence of the two approaches is demonstrated by showing that 
calculated data for Ta II yield the same spacing distribution 
as obtained from observed data for Ta II and Re I combined. 
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These are complex spectra in which a fully developed repulsion 
effect is present. A similar demonstration of equivalence is 
carried out for spcictra of Ru I and Mo I, where the repulsion 
effect is in an intermediate state of development. The 
results also indicate that numbers easily evaluated from the 
radial parameters of the theory will indicate roughly the 
degree of repulsion, replacing to some extent the need for an 
explicit calculation of the spacing distribution. 

Statistical dynamics of simple cubic lattices. Model for the 
study of Brownian motion II, R. J. Rubin, /. Math. Phys. 2, 
No. 3, 373-386 (May- June 1961\. 

New results concerning the statistical dynamics of a heavy 
particle in an n-dimensional {nD) cubic lattice are presented. 
In a well-defined sense, the random motions of a heavy 
particle in a ID lattice and a 3Z) lattice are accurately de- 
scribed by Kramers' equation for a free particle and a har- 
monically bound particle, respectively. A related, but not 
independent, result is that the velocity v{t) and position u{t) 
of a heavy particle in a IZ) lattice and a 3Z) lattice constitute 
two dimensional stationary gaussian Markoff processes. It 
is definitely estabhshed that in the case of a 27) lattice the 

is non-Markoffian. In the 



(v{t)\ 
\u{t)) 



station gaussian process 

course of the analysis, several interesting connections be- 
tween solutions of the discrete lattice equations of motion 
and solutions of the corresponding continuum equation of 
motion (the nD wave equation) are uncovered. 

Partial confounding in fractional replication, W. J. Youden, 
Technometric 3, No. 3, 353-358 [Aug. 1961). 
Fractional factorial designs are based on using some fraction 
of the total possible number of experimental combinations 
of the variables. When resources permit, the experimenter 
may obtain both the direct effect of any variable and the 
effect as modified by the other variables. The latter infor- 
mation may be sacrificed if fewer combinations are exi)lored. 
This paper presents an intermediate possibility without 
increasing the number of combinations. The direct and the 
differential effects may be identified and estimated at some 
sacrifice in the precision of the estimates. 

Exact and approximate distributions for the Wilcoxon sta- 
tistics with ties, S. Y. Lehman, ,/. Am. Stat. Assoc. 56, 293-298 
(June 1961). 

The exact distribution of the Wilcoxon statistic was calculated 
for five specific cases where both samples contained five 
observations, some of which were tied. The five specific 
cases were obtained from an actual experiment. The exact 
distributions for each of these five cases are compared with 
some approximations at points near 1, 5, and 10 percent. 
Critical values and associated probabilities have been cal- 
culated for the exact and approximate distributions. 

Some geometrical theorems for abscissas and weights of 
gauss type, P. J. Davis and P. Rabinowitz, /. Math. Anal. 
Appl. 2, No. 3, 4^8-437 (June 1961). 

Several theorems which relate Gaussian weights and abscissas 
asymptotically are projected and others are conjectured. 
Those conjectured have been established numerically. 
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